Evidence for breakdown of the DGLAP description in diffractive DIS at HERA 
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HERA data on diffractive DIS show deviations from twist 2 DGLAP predictions below 
Q 2 ~ 5 GeV 2 at low pomeron which may reach up to 100 %. These deviations are 
consistent with higher twists effects extracted from the saturation model. It is a first direct 
evidence for the higher twists in DIS. This finding affects determination of the diffractive 
parton densities that are used for the predictions at the LHC. 

I. INTRODUCTION 

O : 
<N ; 

A large fraction of the deep inelastic scattering (DIS) HERA data are produced due to the processes 

I— » 

OO 



of diffractive dissociation of virtual photons y*P ~ > Xp (DDIS). In those reactions, viewed from the 
proton rest frame, the virtual photon y* fluctuates into the strongly interacting debris which scatter 



^ off the proton target and eventually create a final hadronic state X. The proton target remains intact. 

The conventional description of DDIS is based on the leading twist DGLAP evolution equations which 
J|h characterize the QCD hard scale dependence of the diffractive parton distribution functions (DPDFs). 

This approach is justified by Collins factorization theorem [1]. Despite of clear success such description 

^ faces an important limitation that follows from neglecting higher twists. The higher twists contribution 

^ ! becomes relevant below some energy scale which depends on the process. In the case of the inclusive 

, 1 ' DIS, a leading twist description of the data is reasonable down to the photon virtuality Q 2 ~ 1 GeV 2 OO, 
^sD 

O however, in the case of DDIS such description breaks down already at higher scale Q 2 ~ 5 GeV 2 . In this 
<N ; 

^ presentation we would like to show, that the deviations of the extrapolated DGLAP description from the 
• i-h DDIS date at low Q are consistent with emergence of the higher twists contribution. This provides a 

XI ■ 

Sh first evidence of the higher twists effects in DIS data and opens a new window for studying the physics 

a • 

beyond leading twist, both at the experimental and theoretical level. 

II. CROSS SECTION AND THE BREAKDOWN OF TWIST-2 DESCRIPTION 

The DDIS is an quasi-elastic electron - proton scattering process e(k)p(P) — > e(k J )p(P')X(Px) in 
which the final hadronic state X with four-momentum Px is separated in rapidity from the proton, that 
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• ZEUS(LRG) DGLAP — DGLAP + Twist 4 



- DGLAP + Twist 4+6 



DGLAP 

DGLAP + MSS-Sat twist 4+6 
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FIG. 1: Left panel - the % 2 /d.o.f. for NLO DGLAP and NLO DGLAP + HT fits to ZEUS LRG data (3D with 
Q 2 > 2^ in . Right panel - the LRG ZEUS data for £o? (3) at low Q 2 compared to a DGLAP fit OD and the DGLAP 
tit with included twist-4 and twist-4 and 6 corrections from the saturation model. In yellow (gray) — the region 
of |3 where the correction due to qqgg may be neglected. 



scatters elastically. The ^-integrated ep cross-section reads: 

do 



(1) 



d$dQ 2 dt, [3£ 4 

where the invariants read y = (kq)/(kP), Q 2 = -q 2 , \ = (Q 2 +M] ( )/(W 2 + Q 2 ) and t = (P' -P) 2 . 
The quantity W 2 = (P + q) 2 is the invariant mass squared in photon-proton scattering, and M\ is the 
invariant mass of the hadronic state X. The reduced-cross-section may be expressed in terms of the 
diffractive structure functions of (P, Q 2 ,£,) = +F^ 3 \ whereas the structure functions T,L may 
be, respectively, expressed through transversally and longitudinally polarized 7* - proton cross sections 
F?P = (Q 4 /^ 2 a em mdalydM 2 . 

In recent analysis Q3D the ZEUS diffractive data were fitted within NLO DGLAP approximation. A 
satisfactory good description was found only for Q 2 > Q 2 mjn = 5 GeV 2 . The ZEUS fits were performed 
above Q? njn and then extrapolated to lower photon virtualities. The deviations of the fits rapidly grow 
with decreasing ^ and Q 2 reaching 100 percent effect at the minimal Q 2 = 2.5 GeV 2 and \ ~ 4 ■ 10~ 4 . We 
confirmed this result calculated % 2 /d.o.f. for subsets of ZEUS LRG data with Q 2 > Q 2 nin and (3 > 0.035 
Q4J] (see Fig. 1 , left panel). The cut off in P serves to reject contributions from highly resolved projectiles. 
It is clear from this discussion that the leading twist DGLAP evolution is unable to describe the DDIS 
data below Q 2 ~ 5 GeV 2 and at the low 



III. ESTIMATION OF THE HIGHER TWIST CONTRIBUTIONS 



The large energy limit of the DDIS scattering may be described within the framework of the colour 



dipole model [|5i 



]. In this approach the y*P process is factorized into an amplitude of photon fluctuation 
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into the partonic debris and then scattering of these states off the proton by the multiple gluon exchange. 
We take into account the contributions from the fluctuation of the photon into a colour singlet quark- 
antiquark pair qq and into gg-gluon triple. We neglect the contribution from the higher Fock states. 
This gives the /-integrated fp cross section do[j/dM\ = do q ^ T /dM\ + dof^j / ' dM\. Assuming an 
exponential /-dependence of diffractive cross-section, one finds for the qq component 

d<%T 1 f d 2 p f\ Jp 2 

~ 1V1 X j 

f spin 

where bo is a diffractive slope, zz = z(l — z) and the first sum runs over the three light flavours. The 
second sum of © means summation over massless (anti)quark helicities (h)h in the case of longitudinal 
photons whereas for transverse photons there is an additional average over initial photon polarizations 
X. The squared photon wave functions reads, e.g. D7|] 



(2) 



V /- (Q z (Q zf) = ^ *&i*f*o{*)*6&) (L) 



N / (n , ?w/* (n , -z>\ — — L 

spin 

where Kq,K\ are McDonald-Bessel functions and £ = y/ziQ. We use the GBW parametrization [|6|] 
for the dipole-proton cross section o^(r,^) = Oo(l — exp(— r 2 /4R^)). The saturation radius in DDIS 
depends on = (^/x ) X/2 GeV 1 and a = 23.03 mb, X = 0.288, x = 3.04 • 10~ 4 . 

The contribution of the qqg component of y* is calculated at P = and in the soft gluon approxima- 
tion (the longitudinal momentum carried by a gluon is much lower then carried by the qq pair). This 
approximation is valid in the crucial region of M\ ^> Q 2 or |3 <C 1 , where the deviations from DGLAP are 
observed. The correct [3-dependence is then restored using a method described in USD , with kinematically 
accurate calculations of Ref. [9J. With these approximations one obtains: 

dofj 1 AUx,q2 



m X J f spin J0 



dM 2 \6nb D 2% 2 
r r 2 

Nqq g (r 01 ) = / rfVoz^V (Noz + Nn-NmNii -N i) 2 

J r 02 r 12 

where Nij = N(rj — ri), roi,ro2, r\i = ^02 — ^01 denote the relative positions of quark and antiquark (01), 
quark and gluon (02) in the transverse plain. The form of N 2 qg follows from the Good- Walker picture 
of the diffractive dissociation of the photon [10]. The factor 1/M| is a remnant of the phase space 
integration under the soft gluon assumption. The twist decomposition of © is performed using the 
Mellin transform in the roi variable: 

<>°"S N c a 2 a s [ds/AQiy 2 
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The Mellin transform of A^(roi) can be done in two steps. First one defines new integrals 

/ j2„ f,.2 \S~\ I j2„ '01 fAT , AT AT AT ^2 » 

v 01 



h = — / d 2 m (r 2 0i y- 1 [ J 2 r 02 4V \(No2 + N n ~ N 02 N l2 ) 2 - N { 

71 J J r 02 r 12 

f0 2 V f f r 2 

h = / d 2 r 0l (r 2 0i y- 1 / J 2 r 02 ^OJ_2^oi[^o2+M2-^o2^i2-^Voi]. 



(6) 



where N 2 ^ g {s) =I\—I 2 . The integral l\ can be performed exactly, 

/i=Jt(<2o^) 2i 2 1+i (2 1+i -l)r( S )x[^-3 J F 2 (l,l,l-*;2,2;-l>], fl, = £^ 

and for h we use the large daughter dipole approximation ro2 3> roi , ri2 ~ ^02 and obtain, 
/2=^(2o^) 25 2 1+2 'T(5) |i-2 1 - i ' + 3- i '+^ (l +s,l+s;2+s;-^\ 



(7) 



(8) 



The twist decomposition follows from © as a sum over residues at the negative integer values of s. 
Accuracy of this approximation is at the level of 5 per cent. 

IV. DISCUSSION 

In Fig. 1, right panel we compare selected results with data: the extrapolated DGLAP results, 
DGLAP plus twist-4, DGLAP plus twist-4 and twist-6. The saturation model results are obtained using 
the original GBW parameters X and Go, and three massless quark flavours. In our approach we mod- 
ified the GBW parameter xq to ^0 = 2xq in order to account for the difference between Bjorken x and 
^, the variables used in GBW dipole cross-section in DIS and DDIS respectively. We chose a s = 0.4 
that provides a good description of data. The conclusions from the analysis and from the Fig. 1 are the 
following: (i) at twist-2 the DGLAP fit and the twist-2 components of the model are reasonably consis- 
tent, but all fail to describe the data below Q 2 = 5 GeV 2 and at low ^; (ii) a combination of the DGLAP 
fit and twist-4 and twist-6 components of the model gives a good description of the data at low Q 2 . 
Inclusion of the higher twists terms improves the fit quality in the low Q 2 region (see the dashed curve 
at Fig. 1, left panel). Indded, the maximal value of % 2 /d.o.f. ~ 1.5 at Q 2 mjn = 2 GeV 2 is significantly 
lower then % 2 /d.o.f. ~ 3 of the DGLAP fit. However, it is important to stress that a truncation of the 
twist series (up to twist-6) is required to have a good description of the data. Such truncation, however, 
may be motivated in QCD. Let us recall that in BFKL, at the leading logarithmic approximation, only 
one reggeized gluon may couple to a fundamental colour line. Since DGLAP and BFKL approximations 
have the same double logarithmic (lnxln Q 2 ) limit, one concludes that also in DGLAP couplings of more 
than two gluons to a colour dipole is much weaker than in the eikonal picture. Thus one can couple only 
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2 gluons to a colour dipole and up to four gluons to qqg component (two colour dipoles in the large N c 
limit) without violating BFKL constraint. This means that one may expect a suppression beyond twist-8 
if only the qq and qqg components are included in the calculations. This qualitative argument requires, 
obviously, further detailed studies. 
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